A MODEL FOR THE DEGRADATION OF POLYIMIDES DUE TO 

OXIDATION 
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SATISH KARRA AND K. R. RAJAGOPAL 

Abstract. Polyimides, due to their superior mechanical behavior at high temperatures, are used 
in a variety of applications that include aerospace, automobile and electronic packaging industries, 
as matrices for composites, as adhesives etc. In this paper, we extend our previous model in 
[S. Karra, K. R. Rajagopal, Modeling the non-linear viscoelastic response of high temperature 
polyimides, Mechanics of Materials, In press, doi:10.1016/j.mechmat.2010.09.006], to include 
oxidative degradation of these high temperature polyimides. Appropriate forms for the Hclmholtz 
potential and the rate of dissipation are chosen to describe the degradation. The results for a 
specific boundary value problem, using our model compares well with the experimental creep 
data for PMR-15 resin that is aged in air. 
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l. Introduction 

Polyimide and polyimide composites are used in a variety of applications due to the high glass 
transition temperature of above 300°C. These polymers and their composite components undergo 
degradation in a variety of ways including degradation due to oxidation. Thus, there is a need 
to understand how the mechanical behavior of these materials is affected by oxidation. Several 
experimental studies have been carried out which show that there is: (a) wei ght loss in the poly 



imides, and (b) an oxidized l ayer is formed on the surface o f the material (see lBowles et al.l ( 12001 



2003 ). also see references in Pochiraju and Tandon ( 20091 ) ) due to oxidation. The loss of weight 
due to oxidation is observed to be due to chemical bond breakage and escape of volatile lower 
molecular weight gaseous products. In addition, it has been observed that the brittle oxidized 
layer formed on the surface of the polyimide acts as a crack initiation site, which leads to the 
failure of the materials. These cracks als o provide more surf ac e area for further degr a datio n 
and dama g e due to oxidation. Recently, Tandon et al.l ( 2006 ). Pochiraju and Tandon ( 2009 ). 
Roy et al.l ( 120081 ) have looked at oxidative degradation of polymer composites from a modeling 
perspective. However, most of the works either do not consider the coupling between chemical 
reactions and deformation or assume that the coupling is between the small strain in a linearized 
elastic solid model (which does not correctly describe the mechanical behavior of these high 
temperature polymers since it has been experimentally sh own that they exhibit non-linear vis- 
coelastic response, see Falcone and Ruggles-Wrenn ( 20091 )) and an advection-diffusion- reaction 
equation. 

A thermodynamic framework that considers the coupling between chemical reactions (includ- 
ing stoichiometry and chemical kinetics) and deformation of polyimides that show non-linear vis- 
coelastic response, is needed. Such a framework can also be used in modeling similar coupling in 
areas like asphalt mechanics, biomechanics a nd geo-mechanics . Some of the earlier works in areas 
of stoichiometry and thermochemistry are by lPrigogind (119671 ). Ide Ponder and Van Rysselberghe 
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In this paper, we shall extend our cons titutive theory that h a s been used to model the non 



linear response of viscoelastic solids (see Karra and Rajagopal ( 201dh ) to include degradation 



due to chemical reactions (specifically, oxidation). This theory is based on the thermodynamic 
framework of Rajagopal and co-workers (we refer the reader to lRajagopal and Srinivasal (120041 ) 
for details of this framework) that has been shown to be able to capture a plethora of phenomena. 
We extend our previous work by introducing a variable a that represents the extent of oxidation 
in the polyimide. Our approach should not be thought of as merely other internal variable 
theories that are in vogue; we are able to assign a clear meaning to this variable and thus it is 
a variable that goes towards specifying the state of the body . The forms for the Helmholtz ip 
potential and the rate of dissipation £ in lKarra and Rajagopal! (120 lOf ) are modified to incorporate 
the chang es in the response of th e body due to oxidative degradation. Our approach is similar 
to that of Rajagopal et al.l ( 2007 ) who have modeled the degradation due to deformation and 
chain scission in polymers using a variable to quantify the degradation. 

The current paper is organized as follows. The preliminaries that are required are doc- 
umented in section (J2J). In sections (I3.ip . f 13 ■ 2 j) . the constitutive relations for the degrada- 
tion due to oxidation are derived. In section ( 13. 3p . the predictions of our proposed model 
are compared with the experimental creep data for oxidative degradation of PMR-15 given in 
Ruggles-Wrenn and Broeckertl (120091 ). In the final section, we make some remarks concerning 
the limitations of the approach, the scope for its improvement and future work that needs to be 
carried out. 



2. Preliminaries 

Let Kr(B) and n t (B) denote the reference configuration and the current configuration, respec- 
tively. The motion x KR * s defined as the one-one mapping that assigns to each point X G Kr, a 
point x G K t , at a time t, i.e., 

x = XKR (X,t). (2.1) 

The mapping XK R {X,t) is assumed to be sufficiently smooth and invertible. Let K p (t) be the 
configuration instantaneously reached by the body upon removal of the external stimuli. We 
shall call this configuration the natural configuration corresponding to n t (see figure ([T])). Let F 
be gradient of motion x KJS (JC,t) (usually referred to as the deformation gradient), defined by 

f = ^Xkr / 22 n 
dX ' y 1 

and let the left and right Cauchy-Green tensors be defined through 

B = FF T , C = F T F. (2.3) 

The velocity is defined by 



(2.4) 



reference configuration 



current configuration 




'«p(t) 

natural configuration 
Figure 1 . Illustration of the various configurations of the viscoelastic solid. 

Let F Kp[t) be the gradient of the mapping from K p (t) to , and let G be defined by 



(~1 Tp _ jp — 1 El 



(2.5) 



Similar to 02 .3p . we shall define the left Cauchy-Green stretch tensors B G and B p t t \ through 

B G = GG T , B pit) = F Kp(t) F T Kp(ty (2.6) 
We shall also define the velocity gradients 

L = FF\ L G = GG~\ L p = F Rm F~l {t) , (2.7) 
and their symmetric parts by 



Di = - [Li + L i ) , i = p(t), G or no subscript. 



(2.8) 



Also, we define the principal invariants through 



I Bj =tr(B,), II Bj = -{[tr(B,)] 2 -tr(£f)}, = det(B,) l = G,p(t), (2.9) 

where tr(.) is the trace operator for a second order tensor and det(.) is the determinant. Now, 
from (J23D 



F - F Rp{t) G + F Kp(t) G 
FF^ 1 = F K ...GG^F-] , + F K m G 
=4> £ = Lp ( i) + F Kp(t) L G F-i (t) , 



(2.10) 



4n general, it is the mapping from the tangent space at a material point in n p (t) to the tangent space at the same 
material point in n t . 
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where (.) is the material time derivative of the second order tensor. In addition, 

B P (t) = F Kp(t) F T + FF T Kp{t) 



L p{t) B p(t) + B p(t) L p(t)i 



(2.11) 



and similarly 

B G = L G B G + B G L T G . (2.12) 
Hence, from (I2.10P and ( 12. lip , we have 

B pi t) = LB p{t) + B p(t) Lp (t) - F Kp{t) (L G + Lq) (2.13) 

and so 



v 



B p{t) =-2F Kp(t) D G Fi p(t) , (2.14) 

v v 
where (.) is the usual Oldroyd derivative defined through A: = A — LA — AL . The local form of 

the balance of mass, linear momentum, angular momentum (in the absence of internal couples), 

and energy are given by 

g = — gdiv(v), (2.15a) 

gv = div (T T ) + gb, (2.15b) 

T = T T , (2.15c) 

ge = T.L + gr - div(q). (2.15d) 

where T is the Cauchy stress, g is the density, b is the specific body force, e is the specific internal 
energy, r is the radiant heating, q is the heat flux, div(.) stands for the divergence operator in 
the current configuration. 

The local form of the second law of thermodynamics (for non-isothermal processes) is: 

T ■ D — g%l) — gsQ — Qh ' = g9( := £ > 0, (2.16) 

where ip is the specific Helmholtz potential, g is the density, 9 is the temperature, s is the specific 
entropy, q h is the heat flux, ( is the rate of entropy production and £ is the rate of dissipation. 
The kinematics presented in this section are sufficient for the work that follows. 

3. Constitutive assumptions 

3.1. General results. We shall assume that the viscoelastic solid is isotropic and incompressible 
with the specific Helmholtz potential of the form 

tp = vp(B p(t) ,B G ,6,a) = ^(l Bp{tV llB v(tV lBG^ l B G ,0,a), (3.1) 

where a is a variable that accounts for the extent of oxidation, a equal to zero implies the 
body is in its virgin state and a = 1 means that the material is completely oxidized, and no 
further oxidation is possible. The rate of change of the variable a is related to the rate of 
the oxidation reaction that takes place in the polyimide. In general, o: can be taken to be a 
tensor which represents the degree of oxidation in different directions i.e., anisotropic oxidation. 
Although it is seen in experiments that oxidation mostly occurs on the surface of the polyimide, 
we shall assume that oxidation occurs at every point in the body. One can model the motion of 
the surface of oxidation in the polymide using a mixture theory approach such as that used by 
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Rajagopal and Taol fjl995h which can take into effect the diffusion of a singular surface, which in 
the case of the problem under consideration would be the surface that separates the region of the 
virgin and oxidized body; however such an approach would make the problem too complicated 
to be amenable to a meaningful study of a initial-boundary value problem. We are also assuming 
that oxygen, polyimide and the products of oxidation together constitute a constrained mixture 
i.e., there is no relative velocity between these constituents, and hence our approach does not 
capture the diffusi on process. I n order to model the diffusion phenomenon, one can follow the 
approach shown in Karra ( 201dl ). 

Next, assuming that the elastic response from the current configuration n t to the natural 
configuration k p m is isotropic, without loss of generality, we choose k p ^ such that 



F = V 

■* K P (t) r K P (t)' 



(3.2) 



where V K , t . is the right stretch tensor in the polar decomposition of F K , . We shall also assume 
that the total rate of dissipation can be split additively as follows 



rwi T-j ; a c ^ n Qh ; grad(fl) 

T ■ D — gip — gsu = £ m , d > 0, = £ c > 0, 



(3.3) 



where is the rate of dissipation due to the conversion of mechanical working into thermal 
energy and due to degradation, £ c is the rate of dissipation due to heat conduction. Now, if we 
constitutively choose 



g h = -/C(0)grad(0), K{9) > 0, 



(3.4) 



where K, is the thermal conductivity, then (13.31) ^ is automatically satisfied. 
Now, 
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(3.5) 



and using ( T2~TT]) . fP2~T2|) along with in ( 133|) . we obtain that 



if) = 2 
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(3.6) 



Next, we shall assume the rate of dissipation £ m ^ to be of the form 



£ m ,d = €m,d(6, a, oi, Bp(t), D G ). 
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(3.7) 



On substituting (I3.6P into ( I3.3p ^. we arrive at 
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Using (l3~9l)-( |3~Tll in (ET5)l . we obtain 



(T - T p{t) ) ■ D + (T p{t) -T g )-D g -q 



dip . 
da 



a 



= £ m , d (^a,d,-B p(t ),.D G ). 
By virtue of the constraint of incompressibility, we have 

tr(-D) = tr{D p{t) ) = tr(D G ) = 0. 
Since, the right hand side of (13. 12ft does not depend on D, using (I3.13p . 

T = pl + T p{t) , 

where p is the Lagrange multiplier due to the constraint of incompressibilit}Q, with 

dip 

(T p(t) - T G ) ■ D G - g—a = £ mid (9, a, a, B p{t) , D G ), 



which can be re-written as 



dip 



[T - T G ) ■ D G - Q-Q^a = £ m ,d(. e > a > B p(t), D G ), 



using (I3~T3|) and (l3~14p . 



(3.8) 



(3.9) 
(3.10) 
(3.11) 

(3.12) 

(3.13) 
(3.14) 

(3.15) 
(3.16) 



The standard method in continuum mechanics to obt ain constraints appeals t o the n otion that the constraint 
response does not work. It has been shown recently bv lRaiagopal and Srinivasa ( 2005 ) that such an assumption 
is in general incorrect. 
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We shall further assume that £ mjC z can be further additively split as follows: 



£ mid (6>, a, B p{th D G ) = £ m (9, a, B 



P (t), 



D, 



a, a) 



(3.17) 



with each of £ m , C,d being non-negative, so that the second law is automatically satisfied. Noting 
that the first term and second terms on the left hand side of (13. 16[) are the contributions to 
dissipation^] due to mechanical working and degradation, respectively, we shall further assume 
that 



(T - T G ) ■ D G 

dip . 



(3.18a) 
(3.18b) 



Now, we shall maximize the rate of dissipation £ m by varying D G for fixed B p m. That is, we 
maximize the function 



$ := Cm + Ai [U — (T — T G ) ■ D G ] + A 2 (I ■ D G ), 
where Ai, A2 are the Lagrange multipliers. By setting, d§/dD G = 0, we get 



(3.19) 



T G + ^I 



Ai 



dD 



G 



(3.20) 



We need to determine the Lagrange multipliers. On substituting ( I3.20p into (13. 16ft . we get 



Ai + 1 
Ai 



9gm / 1 : 

9Dg ' G 



and so (13. 20ft with (13. lip becomes 
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dip dip 
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(3.21) 



(3.22) 



where A := y- is the Lagrange multiplier due to the constraint of incompressibility. 
Finally, the constitutive relations are given by 
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(3.23a) 

(3.23b) 

(3.23c) 
(3.23d) 



The two equations (I3.23ap . (13.23bj) are to be equated and simplified to get the evolution equation 
for B K . This will be shown in the next subsection when we choose specific forms for ip and £. 



3 The term dissipation is used to refer to the mechanical working being converted into energy in thermal form, 
and associated with this dissipation we have entropy production. We shall abuse the use of the term dissipation 
and refer to other entropy producing mechanism such as degradation as also dissipation. 
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3.2. Specific case. We choose the specific Helmholtz potential as 



in 



A s + (B s + c s 2 ) (, 
/x p (l + 7 (0)a) 



0. 



//g (1 + /3(fl)g) 
2^ 



(W " 3) 



+ 



2g 



(I^w " 3 ) + 



(3.24) 



where ftcftp are elastic parameters, S is a reference temperature for the viscoelastic solid and 
the rates of mechanical dissipation, and the dissipation due to degradation as 

U = 77(1 + 5{9)a) (D G ■ B p{t) D G ) , (3.25a) 



D(\\a\ 



(3.25b) 



a 



where rj is the viscosity, |.|| stands for absolute value. Here, /3, 7 and 5 are material parameters 
that depend on temperature. Also, note from fl3.25|) that £ m , ^ are non-negative provided 77, 5, 
D are also non-negative. 
Now, 



s = — 



dt/j 

do 



{B s + 4) + c{{6- 9 S 



din — 



H G a dfi 
~ 1 2o"d9 



(Ifl e " 3) 



\x p a dj 



(I 



2q d6"~ Bp(t) 
The internal energy e is given by 

e = ljjJ r 0S 



6) 00- 



(3.26) 



A s - B S 9 S + 4 
/7 



0, 



|(^ 2 -^ 2 ) + |f(lB G -3) 



1 + a /3 



3) 



1+a 7- 9 



and the specific heat capacity C v is 

de~ 1+2 2 Q 



89 



F-e 



OF 
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Also, fl3T23ij) . f l3~23b|) reduce to 

T = pi + ft p B p( t), 

T — XI + ft G B G - 



(B p ( t )D G + DaB p {t)j 



(3.27) 
(3.28) 

(3.29) 

(3.30a) 
(3.30b) 



where /x p = fi p (1 + /3(6 l )a), /i^ = ftG (1 + 7(0)cO, 77 = 77 (1 + 5(^)a). We also note that we chose 
the functions for the material moduli (ft p , ft G , 77) such t hat they increase as a goes from t o 
1. This is consistent with the experiments (see figure 5 in Ruggles-Wrenn and Broeckert ( 20091 )) 
where it seen that the elastic modulus increases with aging . We further note that such a choice 
of functions for the material moduli is different from what iRajagopal et al.l ( 120071 ) have chosen 
in their work. 



From (ETSQjl 



77 

(p - A)/ + ppBpft) = Ji G B G + - (B p{t) D G + D G B p{ t)) 



and so by pre-multiplying the above equation by S ,^ and taking the trace, we get 

Jkfa(B~{ t) B G ) - 3/Zp 



(p-A) 



Using (I3.32p in (I3.3ip . we arrive at the following equation: 



ft G U(B^ t) B 



3wj 



which can be re-written as 
^^(B-^g) - 3/2 p 

. tr(B; ( J } ) 



/ 

I + fipBpty = p G B G + - (B p ( t )D G + D G B. p (t)) , 



(3.31; 



(3.32) 



(3.33) 



J + // p Bp( t ) = ii G B G 



-?[V 



p(t) *>p(t) V « p(t) I" V Kp(t) £»p { t) V p{ t) I , 



where we have used (12. 14)) and (13.331) . Now, ( I3.23dl) reduces to 
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We shall assume that 



F(a) 



k{9) 



-a, 



where k is a non-negative constant, then (I3.35p . for n — 1 reduces to 



(I 



2 v " Sg ~ y ' 2 " d(l-a) 

Notice that the first two terms on the left hand side of (13.371) represent the dependence of the 
extent of oxidation on the deformation of the material. 

Thus, with the current choice of the specific Helmholtz potential and the rate of dissipation, 
we arrive at the following constitutive equations: 

T = pI + n p (l + P(9)a)B p(t) , (3.38) 

where the evolution of the natural configuration is given by 

(1 + l(0)a) tr(B^)Bo) -3^(1 + 



^(V,"3)-* 



D\\a\ 



(3.34) 



(3.35) 



(3.36) 



(3.37) 



a 



1 + H P {1 + 0(9)a) B p{t) 



ti (l + -y(9)a) B G 



V p {t) B m 



Vnl {t) + V^ 1 B p{t) V p{t) 



(3.39) 



and the evolution of a is given b y f|3.37p. These constitutive relations reduce to the non-linear 
viscoelastic solid model derived in iKarra and Rajagopall ( 120101 ) when there is no degradation. In 
a general initial-boundary value problem, one has to solve the coupled equations ( 12.15b|) . (13.371) 
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along with f !3.38|) . ( I3.39p . subject to appropriate initial and boundary conditions. In problems 
where temperature gradients are important one needs to also consider the balance of energy 
(l215d]) . 

3.3. Comparison with experimental data. In order to compare the predictions of our model 
with experimental data, we shall consider the problem of uniaxial extension, given by 

1 1 

y 



X 



X(t)X, 



within the context of this model. The velocity gradient is given by 

L = diag < — — r, r 

& |A 2A 2A 

We shall assume that the stretch B p ^ is given by 
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Straight forward calculations using ( I3.33|) give 
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2 A fj B fj 
which can be re-written in the following form: 
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(3.40) 



(3.41) 



(3.42) 



(3.43) 



(3.44) 



Using ( I3.42p in ( I3.38p . and also the fact that lateral surfaces are traction free, it is easy to see 
that 

fl p (b - J , (3-45) 



T 



ii 



and so 



In addition, (I3.37P becomes 



Tn = fi p 1 



2^fB 
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(3.46) 



(3.47) 



We shall also use logarithmic strain (or true strain) e = InA as our strain me asure in what follows. 



We shall compare our model with the experimental data for PMR-15 from lRuggles-Wrenn and Broeckert 
(120091 ). With the given T n and material parameters, (13.461) . (13.441) were first solved using the 



T 

initial condition that 5(0) = A(0) = 1 for a time of -A Then, f[3T44|) . ( 13T45|) were solved till 

^ii 

the end of loading. Since in the experiments the aging in air was done without any load being 
applied, ( I3.37P was also solved without the first two terms on the left hand side using a(0) = 
as the initial condition. The ODEs were solved in MATLAB using the ode45 solver. The follow- 
ing parameters were used for comparing the results predicted by our model to the experimental 
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creep data for PMR-15 (under 10 MPa loading) that has been aged in air for various amounts 
of time (also see figure ©): |i p = 2x 10 9 MPa, \i G = 3.8 x 10 8 MPa, 77 = 45 x 10 12 MPa.s, 
I = 1.2 x 10~ 6 s _1 , (3 = 10, 7 = 0.3, 5 = 0.5. Also, since the experiments were done under 
isothermal conditions, the balance of energy f!2.15dj) was not considered. 



4. Concluding remarks 

A model for the degradation of polyimides due to oxidation has been developed in this paper. 
Our model also accounts for the effect of deformation on the aging due to oxidation. However, 
there is no experimental data to corroborate this part of our model. 

The limitations of our model are as follows: 

(1) Our model cannot predict the diffusion of oxygen and hence one cannot estimate the 
thickness of the oxygen layer on the surface of the resin. 

(2) The weight lost due to oxidation cannot be estimated using our model. For this one 
needs to understand the chemical kinetics. Once the chemical kinetics are established, an 
approach similar to ours can be used to couple these reaction kinetics to the deformation 
of the polymer. Our work in this paper can be viewed as a first step towards this end. 

Directions for future work follows: 

(1) Further experiments need to be carried out to get the creep and stress relaxation of 
different polyimides under different loading conditions and for various amounts of aging 
in oxygen. Also, one needs to perform experiments where simultaneous deformation and 
aging takes place to see how they are coupled. 

(2) One can easily extend our work to include anisotropy in oxidation by introducing a tensor 
(a), and also anisotropic response of the body by modifying the Helmholtz potential. In 
addition, for modeling oxidative degradation in polyimide composites, one can also use 
different variables a m , atf that represent the degradation du e to oxidation in t he po lymer 



matrix and fiber, respectively. This approach is similar to iBaek and Pence! (120091 ) who 
have used different variables to represent the degradation due to swelling in the composite 
matrix and fiber. 

(3) Our framework in this paper can also be used to model oxidation or other chemical 
reactions in materials like asphalt / asphalt derivatives, and biomaterials. For this, one has 
to accordingly modify the terms in the Helmholtz potential and the rate of dissipation 
due to deformation. 
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